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An overlooked conservation law for near-inertial waves propagating in a steady back¬ 
ground flow provides a new perspective on the concentration of these waves in regions 
of anticyclonic vorticity. The conservation law implies that this concentration is a direct 
consequence of the decrease in spatial scales experienced by an initially homogeneous 
wave field. Scaling arguments and numerical simulations of a reduced-gravity model of 
mixed-layer near-inertial waves confirm this interpretation and elucidate the influence of 
the strength of the background flow relative to the dispersion. 


1. Introduction 

Near-inertial waves (NIWs) are ubiquitous in the o cean. They con t ribute strongly to 


surface mixing, and hence to biological activity, e.g. iGranata et al\ (jl995l ). They also 
propagate to depth where they eventually di ssipate, thus participatin g in deep vertical 
mixing and the global overturning circulation ( Ferrari &; Wunschlll2009l ). The propagation 


of NIWs in heterogeneous flows has motivated a great deal of wor k. One o f the main 
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Danioux et al . I 2008f ) and observational data ( Knnze fc Sanfordlll984t Elipot et al 


.e., m 


regions of negative relative vorticity in the Northern hemisphere and positive relative 
vorticity in the Southern hemisphere. 

Various explanatio ns have been advanced for this phenomenon. Using a WKB ap¬ 
proach, iKunzel (Il985t ) showed that refraction by a background flow leads to an ‘effective’ 
inertial frequency f e shifted from the local Coriolis frequency / by C/2, where C is the 
local relative vorticity: f e = f + C/2. Because the range of allowed wave frequencies uj 
satisfies |/ e | < c j < TV, where N is the stratification frequency (usually much larger than 
|/|), it is larger in regions of negative vorticity in the Northern hemisphere. Therefore, 
anticyclonic regions are expected to be more energetic because (i) NIWs propagating 
poleward can enter regions of negative vorticity regions, but not regions of positi ve vor¬ 
ticity ; and (ii) frequencies lower than / present in the forcing can become resonant ( Kunzel 
Il985f b This explanation is subject to caution, however, because the initial length scales 
of NIWs, set by atmospheric storms, are typically much larger than the length scales of 
geostrophic flows (the ratio of tropospheric to oceanic radii of deformation is about 10), 
thus invalidating the assumption underlyin g the WKB approximation. 

Alte rnative explanations were offered bv lYoung fe Ben Jelloull (1997) an d Klein et~al\ 
(|2004h . These rely on the NIW model developed by Young fe Ben Jelloull (jl997h . here¬ 
after refer to as YBJ model, which exploits the small frequency spread of NIWs near / 
and makes no assumption about their spatial scales. These explanations rely on strong 
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assumptions: strong dispersion for lYoung fe Ben Jelloull ’s. short-time and restrictions on 
the spectrum of the vorticity field for lKlein et a/J ’s. In this note, we revisit the issue and, 
taking the YBJ model as a starting point, show that a so-far overlooked conservation 
law provides a robust argument for the concentration of NIWs in anticyclones. Section [2] 
derives the conservation law in the simple case of a reduced-gravity shallow-water model, 
the YBJ approximation of which is obtained in Appendix \K\ (The extension to a contin¬ 
uously stratified fluid is straightforward when the background flow is assumed barotropic 
so that NIWs can be expanded in vertical modes.) Section [3] demonstrates analytically 
and numerically how NIW concentration in anticyclones stems from the conservation 
law. Some conclusions and perspectives are offered in section 0J 


2. Model and conservation laws 

2.1. YBJ model and analogy with the Schrodinger equation 

We study the propagation of NI Ws in a steady geostrophic flow with a reduced-gravity 
shallow-water model (e.g.. ICushman-Roisinl 1994 ). This slab-model can be thought of as 
representing the dynamics of NIWs confined in a mixed layer capping an abyssal layer 
where the only motion is the imposed geostrophic flow, assumed to be identical in both 
layers. For NIWs, the horizontal velocity ( u , v ) can be written in terms of a complex 
amplitude M according to u + iv = Me~ 1 ^ 1 . The slow (compared to / -1 ) time evolution 
of M is governed by the YBJ equation 


dtM + M) — i^AM ■ 


AY 

-iJfM = 0 , 


( 2 . 1 ) 


where ip and Aip are the streamfunction and vorticity of the steady geostrophic flow, 
J is the horizontal Jacobian, and h = g'H/f is a dispersion parameter, with g' and H 
the reduced gravity and average depth of the mixed layer. We assume that h > 0, as in 
the Northern hemisphere where / > 0. The respective terms quantify the effects due to 
advection, dispersion, and refraction. A concise derivation of HUD is given in Appendix 

m 

Two simple facts help in understanding the dynamics of (j2.1D . First, for any constant 
a > 0 the YBJ equation is invariant under the scaling transformation ip aip, h ah, 
t i —y t/a. This makes obvious that the intrinsic dynamics of (j2.1|) depends only on the 
single non-dimensional parameter h/T, say, where T is the amplitude scale of ip. 

Second, without the advection term J(ip,M) the YBJ equation is identical to the 
Schrodinger equation that governs the complex wave function (p(x, y, t) for a single par¬ 
ticle with unit mass and external potential V(x,y): 


d t (p 


. ft. .v 

i—Acp + 1- 


= 0 . 


( 2 . 2 ) 


2 T ' ' ” ft 

Here h is Planck’s constant divided by 27r. Comparing and (IT2l) and identifying h 
with h shows that the effective potential in (|2.1[) is 

Aip> 


V = h- 


(2.3) 


Clearly, regions of higher V repel the particle whereas regions of lower V attract it. Hence, 
if the advection term is negligible, then the mathematical analogy between the particle 
probability density \<p\ 2 and the inertial wave kinetic energy density | M| 2 readi l y im¬ 
plies that cyclones repel inertial waves whilst anticyclones attract them ( Balmforth et all 
Il998h . For the special case of an axisymmetri c anticylonic vortex, t his is confirmed by 
the existence of axisymmetric trapped modes ( Llewellyn Smithlll999l ). 
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Of course, if the advection term cannot be neglected then the simple analogy with 
the Schrodinger equation breaks down. Hence, the main task is to understand how the 
advection term alters the basic Schrodinger dynamics as a function of h/^f. For this it 
becomes crucial to study the full set of conservation laws associated with (|2.ip . as we 
shall do now. 


2.2. Conservation laws 


Multiplying (j2.ip by M* and adding its complex conjugate gives 

V • F = 0, 


d t ^\M\ 2 + J(^\M\ 2 ) 


(2.4) 


where F = M* VM)/4. In a finite domain with suitable boundary conditions 

(periodicity o r M = 0), integrating (12.41) gives the conservation of NlW-kinetic energy, 
as derived by lYoung fc Ben Jelloull (jl997h , 

Jt //\\ M \ 2(ixd y = °- ( 2 - 5 ) 

There is another conservation law associated with (|2.ip and not mentioned in Young fe Ben Jelloull 
(1997). It is derived by forming the combination M* (12.11) — M t { 2. ID * and integrating the 
result over the domain. Using properties of the Jacobian and integrating by parts, this 
gives 


J t ^+h + h) = 0, 


(2.6) 


where 


h = 


Jj ihipj{M*,M)dxdy, h = JJ y|VM| 2 ^dy, I 3 = JJ 


h^J-\M\ 2 dxdy. 


The terms I \, I 2 and Is stem directly from the advection, dispersion and refraction terms 
in the YBJ equation. The overall factors of h are included for two practical reasons: to 
keep the values of the invariants comparable when h is varied, and to highlight the 
appearance of the effective potential V from (E2) in J 3 . The consequences of the new 
conservation law d!zi)ii are discussed in more detail in the next section. In the absence of 
a geostrophic flow, ii, Is are identically zero, and I 2 can be recognised as the scaled NIW 
potential energy averaged over the fast time scale. This is not unexpected, as explained 
in Appendix [Bj In the presence of a steady flow, I\ + I 2 + Is can be interpreted as an 
energy in that its conservation is associated with the time invariance of (EU. Unlike 
(EU), h + I 2 + ^3 is not conserved for arbitrary time-dependent flows. However, when 
flow and NIWs evolve in a dynamically consistent manner, an analog ous conservation 
law holds that accounts for energy transfers between flow and NIWs ( Xie fe Vannestel 
2015lh 

We remark that, in a steady geostrophic flow, differentiating 42U with respect to time 
shows that M t satisfies the same equation as M, and hence the same conservation laws. 
In particular, 


ill^ 


l dxdy = 0 


(2.7) 


means that the root-mean-square magnitude of M t is constant. 


3. NlW-concentration in anticyclones 

Non-dimensionalizing EID using x = Lx ', ^ h = Th/, and t = (I/ 2 /T)t', 

with L the typical length scale of the geostrophic flow, gives an identical equation for 
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the primed variables. Again, this makes obvious that hi = h/T is the only relevant 
parameter; in £13.31 we conduct simulations with different values of h/T. This parameter 
is the reduced-gravity shallow- water equivalent to the param eter T = ^t/h = 1/h' used 


Young fc Ben Jelloull ( 19971 ) and iBalmforth e£~a/1 ( 1998 ). on which they base their 


‘strong dispersion’ (h/T 1) and ‘strong trapping’ (or ‘strong advection’, h/T <C 1) 
approximations. In the ocean, h/T is highly variable because of varying kinetic-energy 
levels and stratification. For instance, typical values for the North Atlantic might be 
/ = 10 -4 s -1 , g' = 2• 10 -3 ms -2 , L = 50 km; taking H and U in the ranges H G [50,100] 
m, U G [0.005, 0.1] ms -1 gives h/T G [0.2,8]. 

We will consider a specific initial-value problem in which M(x,y,0) = 1. Without 
loss of generality this represents an eastward NlW-momentum deposition by a storm: 
because storm scales are typically much larger than ocean eddy scales, a homogeneous 
initial condition for M is appropriate. 


The conservation law involves Is , which is proportional to the covariance between 

\M\ 2 and A-0; it is therefore relevant to the concentration of NIWs in anticyclones, which 
corresponds to Is < 0. At t = 0, I\ = h = Is = 0 (assuming no net vorticity - a given 
with periodic boundary conditions). The development of spatial heterogeneities in the 
M- field must lead to an increase in the positive definite / 2 , which is then compensated 
by h + Is < 0. Of course, if Ii is negligible then Is must take negative values, but in the 
general case it is less clear whether Is behaves in this way. We next provide an asymptotic 
argument that Is becomes negative for all values h/Sfr at short times. This is followed by 
a long-time scaling argument for m as a function of h/T, which predicts that Is < 0 
for large and intermediate values of h/T but not for small values of this parameter. These 
predictions are then checked against numerical simulations in § 13.31 


3.1. Short-time solution 


With homogeneous initial conditions, the first physical effect on NIWs propagation is 
due to refr action, and th e shor t-time behaviour is M(x,y,t) = exp(— itA'i/j(x, y)/2) = 
M{x,y,t) ( Danioux et a/.ll2008 ). A cautious definition for short time here is t <C t s = 
1/max{T/(2L 2 ), h/(2L 2 )}. For such times, the solution can be sought as the expansion 
M = M + M', with M'(x, y, t = 0) = 0 and \M'\ <C 1. Introducing this into (J2.ip gives 

d t M' + J(V>, M') - i^AM' + M) + AM. (3.1) 


Because of the form of M, the short-time behaviour of the right-hand side of m 
behaves as 0(£), hence forcing M'(t ) = 0{t 2 ) for small t. Keeping this in mind, we 
now compare the relative size of terms I\ and I 2 in (j2.6|) . Firstly, because J(M*, M) = 0, 
I\ is dominated by terms of the form ih^ J(M*, M'), resulting in a 0(t 3 ) dependence. 
Secondly, I 2 is dominated by h 2 |VM| 2 /2, yielding a 0(t 2 ) dependence. Hence, for times 
short enough, I 2 I\. NlW-energy concentration in regions of negative vorticity can 
also be deduced from RIM . Injecting the short-time solution M(x, 7 /, t) into F, one finds 
that the amplitude of M obeys 


dt\\M\ 2 



(3.2) 


at short times, consistent with iKlein et~al. (2004)’s short-time solution. Because the 
vorticity field and its Laplacian are anticorrelated, (IQl gives an increase of NlW-energy 
in anticyclonic regions for all values of h/T. 
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3.2. Long-time scaling arguments 

For long times the spatial scales of M need not be equal to those of ip anymore and hence 
the dominant balance between the various terms in (EUi may shift accordingly. However, 
the conservation law (12.71) implies 

Jj \\M t \*dxdy = JJ = 0)| 2 dxdy = JJ (A^/2 fdxdy (3.3) 

and together with (|2.5|) this means that for all times the root-mean-square magnitudes of 
M and M t scale with unity and T/L 2 , respectively. This allows a simple scaling analysis 
of the four terms in dm which after rearranging yields 


L hL_ 

l' ¥72 


(3.4) 


Here l is the long-time spatial scale of M such that VM = 0(1//). We use (|3.4|) to 
determine how l/L may depend on h/T. First, in the ‘strong dispersion’ regime h/T 
1 the only possible balance in (|3.4[) is l/L = ^/h/T, which balances dispersion and 
refraction whilst advection is negligible. As expected, this reduces the YBJ dynamics to 
that of the Schrodinger equation, so I\ is negligible and Is ~ —h- In this regime the 
spatial scale of M is larger than that of the background flow ^ by a factor of yh/T. 

Second, in the opposite regime h/^t <C 1 one scaling possibility is l = L, which balances 
advection and refraction whilst dispersion becomes negligible. However, from (|2.1[) this 
would correspond to an advective dynamics along streamlines in which \M\ is conserved 
whilst its phase continues to change by the refraction as in the short-term solution derived 
in §o This would inevitably lead to the generation of evermore smaller spatial scales 
in M and hence defeat the assumption l = L. We must therefore look at the alternative, 
a balance between the advective and the dispersive terms based on l/L = h/T. The 
advective and dispersive terms then provide a new leading-order dominant balance of 
size T/h 1 to the long-term evolution of the YBJ equation. Consequently, in this 
scenario the refraction is weak and Is becomes negligible and so I\ ~ —I?,. Moreover, the 
spatial scale of M is much smaller than that of -0, by a factor of h/T. 

Note that this small horizontal scale may invalidate the assumption of small h/(fl 2 ) 
that underpins the near-inertial YBJ model (see Appendix E|- Using / = AL/T, we 
rewrite this assumption as T/(/L 2 ) <C h/T in terms of the Rossby number on the left- 
hand side. Hence for small but fixed h/T the near-inertial approximation holds provided 
the Rossby number is small enough. 

Finally, in the intermediate regime h/T = 1 the scaling l/L = 1 makes all terms equally 
important. It stands to reason that in this intermediate regime Is will take moderate 
negative values, consistent with a monotonic transition between its value Is ~ 0 for 
h/T 1 and its negative value Is ~ —12 for h/T 1. However, only in this intermediate 
regime does the length scale of M equal that of ip, which arguably is the best situation for 
effective concentration of NIW energy in anticyclones. This suggests that the intermediate 
regime might be the most effective for this purpose. Of course, this simple scaling analysis 
can only provide a heuristic guide to the full NIW energy dynamics, not least because 
the spatial scales of |M| 2 are not related in a trivial way to those of M. Still, we will see 
that the present scaling arguments do indeed provide a useful guide for understanding 
the numerical simulations. 


3.3. Numerical simulations 

Eq. (|2.1[) is solved numerically on a doubly periodic 256 x 256 grid using a pseudo- 
spectral time-split Euler scheme. A weak biharmonic dissipation is added for numerical 
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stability. The streamfunction 'ip is taken as a single realization of a homogeneous isotropic 
Gaussian random process, with Gaussian correlation function and a correlation length L 
of about 1/5 of the domain size. The chosen vorticity field is shown at the top of figure 
[TJ Because its correlation scale is much smaller than the size of the domain, the results 
presented here are generic. We run three simulations with the same background flow but 
different values of h such that h/T = 0 . 2 , 1 and 10 , representative of the strong advection, 
intermediate and strong dispersion regimes. Simulations are stopped when the various 
terms in (12.61) no longer evolve significantly, which happens around t ~ 0.3 L 2 /h (figure 
[1]). For the intermediate case, using the values given at the beginning of this section, we 
find a time-scale between 4 and 8 days, which is realistic. 

Snapshots of \M\ are shown in figure |T] for very short, short and long times for the 
three values of /i/4L At very short times, for all values of h, \M\ is clearly anti-correlated 
with the vorticity field (or correlated with its Laplacian), as predicted by (|3.2|) . As time 
increases, the variance of \M\ increases, substantially for small h/T but much less for large 
h/T. This is consistent with the scaling \M — 1| = 0(T/h) that holds for short time and 
arbitrary h/T (as follows from integrating (j3.2|) for t < L 2 /h) and for all time and large 
h/T (as follows from a perturbative treatment of (| 2 . 1 |) ). At the end of the simulation, 
M has larger scales than at the very first i nstants for h/T = 10, in accordance with 
the scaling arguments in § 13.21 and also with lKlein et al. ( 2004 )’s ‘truncated Laplacian’ 


solut ion. The spectrum of M does not change much afterwards (not shown), in contrast 
with Klein et al. ( 20041) ’s theory which predicts a continued cascade towards large scales. 
As expected, the evolution is very different for h/T = 0 . 2 . The short-time solution (figure 
ffli) is consistent with a passive-scalar type scenario: after the initial generation of spatial 
scales in \M\ (figure QJi) , the |M|-field is stretched and folded, while its amplitude grows 
due to dispersion. Later on (figured^), NlW-scales much smaller than the flow scales 
develop, consistent with ^3.21 

The behaviour of \M\ is reflected in the evolution of ii, I2 and Is displayed in figured 
The generation of small scales by refraction subsequently modulated by dispersion leads 
to an increase of I 2 - This is balanced by the decrease of Is at short times for all values of 
h/T and for all times for h/T = 10 , as predicted by our scaling argument. Conversely, 
for h/T = 0.2 we obtain I\ ~ —I 2 , indicative of the balance between advection and 
dispersion noted above (in this case, the amplitude of these quantities decrease slightly 
in time because of numerical dissipation). For h/^t = 1, all terms ii, I 2 and Is have 
similar final amplitudes. Note that the amplitude of I2 is roughly the same in the three 
simulations. From the scaling \M — 1| = 0(T/h) mentioned above, we infer that I 2 is 
proportional to T 2 /L 2 , and similarly for Is . This justifies a posteriori the inclusion of 
the overall factor h in the definition of ii, I2 and Is- 

We now turn to the concentration of NIW energy in anticyclones. This is best quantified 
by / 3 /h, that is, the covariance between AtJj/2 and |M| 2 . An alternative diagnostic is the 
correlation between |M | 2 and A?// 2 , namely 


C = 


h 


h (°m If (AW2 ) 2 dxdy) 


1/2 ’ 


where cf 2 m = JJ (|M | 4 — 1 )dxdy (3.5) 


is the variance of the NIW energy. We emphasise that C focuses on the match between the 
spatial patterns of NIW energy and vorticity irrespective of their amplitudes; the more 
intuitive colocation of high values of NIW energy with regions of anticyclonic vorticity is 
measured by Is/h. The evolution of both Is/h and C is shown for the three values of h/T 
in figure [3j At the end of the simulations, the covariance Is/h is largest for intermediate 
values of h/T. Small values of h/T, lead to a strongly heterogeneous \M\ 2 (large cr^, as 
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Figure 1. Streamfunction (a 5 ) and associated normalized vorticity field (b’) of the background 
flow, and evolution of the NlW-amplitude \M\ for h/T — 0.2 (left column), h/T — 1 (middle 
column) and h/T — 10 (right column). Panels (a), (b), (c) correspond to very short times after 
the start of the simulations (1/80 th of the length of the simulations), panels (d), (e), (f) to short 
times (1/ 16 th of the length of the simulations) and bottom panels to the end of the simulations. 


seen on figure Wg), but to relatively weak correlation of |M| 2 and Aip (figure USE dashed 
line); conversely, large values of h/^t lead to a strong correlation (figure [3b, solid line) 
but weak heterogeneity of |M| 2 . The intermediate case displays both a relatively strong 
<j\j (figure Hh) and a strong correlation (figure [3b, dotted line), giving the maximum 
covariance. 


4. Discussion 

The results described in this paper in the context of the reduced-gravity shallow-water 
system demonstrate that the development of spatial heterogeneity in a homogeneous field 
of NIWs due to the interaction with the barotropic vorticity field is inevitably accompa¬ 
nied by a concentration of NlW-energy in anticyclones. We emphasise that this result does 
not make assumptions abo ut the relative importance o f dis persion (such as the str ong- 
dispersion approximation of I Young fe Ben Jelloull ( 19971 ) and lBalmforth et~cd . I ( 1998 )) or 
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h/$> = 1 



h/V = 10 



ht/L 2 


Figure 2. Evolution of the integrals 1 1 (dotted line), I 2 (dashed line) and I 3 (thin solid line) 
in (12.61) for h /T = 0.2 (a), 1 (b) and 10 (c). ii, I 2 and I 3 are scaled by the area of the domain. 
The inset in (a) represents a zoom on very short times. 



(a) 

-0.1 
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-0.3 
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-— h/$ = 0.2 
hj'l> = 1 
-= 10 


0 0.1 0.2 0.3 

ht/L 2 



Figure 3. (a) Spatial covariance 1 3 /h and (b) correlation C between \M \ 2 and A^/2 for 

h/V = 0.2, 1 and 10. 


the nature of the vorticity field and smallness of advection ( Klein et all 2004). Instead , 
it aris es as a consequence of a conservation law associated with I Young fc Ben Jelloull 
( 1997 Vs model. Therefore, the concentration of NlW-energy in anticyclones is a much 
more robust phenomenon than previously thought. In particular, it is the strongest for 
the intermediate regime h/^f = 0(1), when refraction, dispersion and advection are all 
significant. 

We note that (|3.3p leads to an interesting conclusion for h/^t 1. In this limit, M 1 
as explained above so that the refraction term in (EH) becomes approximately iAip/2. 
By (|3.3 p the time-derivative term has the same amplitude as the refraction term in 
(2 .11). in contrast with previous treatments of the s trong-dispersion limit which neglect it 
( Young fc Ben Jellonll997l : lBalmforth et a l 1998). In this case, a more complete solution 
of (|2.ip at short times is 


M = 1 + 




(l-e 


-i/i|k| 2 i/2 


)dk + 0((«v7i) 


(4.1) 


where ^(k) is the Fourier transform of the streamfunction ib at wave number k. From 
631, the adjustment of the initial condition M = 1 towards the ‘balanced state’ with 
M = 1 + i/j/h + is accompa nied by the emission of Q(^f /h )- amplitude waves. 

Altho ugh the latter part is missing from Young fc Ben Jelloull (1997) and lBalmforth et al. 


(Il 998lh 63]) shows the correlation between \M\ and ip (hence the anti-correlation between 
\M\ and the vorticity field) remains true on average. Note th at a slow modulation of (|4.1D 
should be added to describe the long-time behaviour of M ( Young fc Ben Jelloull 1997F . 

In this paper, we make the strong assumption of a steady background flow. The ad- 
vective time scale L 2 /\k typical of the flow evolution can be compared with the 0(L 2 /h) 
time taken for A, I 2 and Is to reach saturation (see figure [2j) to conclude that the as¬ 
sumption of steadiness can be relaxed when h T. When h/^f = 1, saturation occurs for 
ht/L 2 ~ 0.3 so the impact of unsteadiness can be expected to remain weak. In the strong 
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advection regime h/T <C 1 , however, the time dependence of ip cannot be neglected and 
only the short-time solution described in 13. II remains strictly valid. 

Last, the assumption that NIWs initially have much larger horizontal scales than the 
geostrophic flow is crucial for the results reported here. Although this holds in many parts 
of the ocean, NIWs can also be generated at scales similar to those of the geostrophic flow 
(e.g., by a moving hurricane). The study of the propagation of NIWs in a geostrophic 
flow with similar scales is the subject of a forthcoming paper. 

Acknowledgements. This research is funded by the UK Natural Environment Research 
Council (grant NE/J 022012 / 1 ). 


Appendix A. Reduced-gravity shallow-water YBJ model 

We begin by considering the reduced-gravity shallo w-water system line arized about a 
barotropic geostrophic flow (£/, V) = (— r ip y ,'ip x ) (e.g.. iKlein et al\ (|20Q4h ). That is, we 
assume u, v <C C7, V and 77 <C H, where H is the horizontally averaged depth of the top 
layer. We emphasize that, because the geostrophic flow is barotropic (i.e. the same in the 
top layer and beneath), there is no associated interface slope. Under these assumptions, 
the NlW-velocity and layer-depth perturbations obey 


d t u + Ud x u + VdyU + ud x U + vd y U — fv = — g' d x r\ , (A 1) 

d t v + Ud x v + Vd y v + ud x V + vd y V + fu = -g f d y rj, (A 2) 

d t r] + Ud x r] + Vd y r] + H(d x u + d y v ) = 0. (A3) 

Nondimensionalising using (x,y) = L(x',y f ), ip = Tt//, (u,v) = U w (u' ,7/), 77 = HU w /(fL) 7/ 
and t = t' //, and introducing a slow time-scale r = et, we obtain the following equation 
for the complex velocity U = u' + iv' by forming (TOl +i lTOl : 


d t U + iU = —e 


d T U + J{iPM) + 2 %. 


+*—w- 


2 iV> € . € .W 


(A 4) 


where ^ = x + iy and A is the horizontal Laplacian. Primes in (|A4|) have been omitted 
for simplicity. We have assumed e = T/(/L 2 ) ^ g’H/{f 2 L 2 ) < 1; this corresponds to 
assuming a low Rossby number for the background flow and a small Burger number for 
the waves (i.e. waves oscillating at a frequency close to /). The nondimensional version 
of (1A3 [) is 

dty + e J ( 7 ^, T}) + (W^ + ) = 0. (A 5) 

An approximate solution can be sought by expanding U in powers of e: 

U =U (0) + eU w + 0(e 2 ). (A 6) 

The leading order solution can be written as 

Z/ 0) = M(x,y,T)e- u , (A 7) 


where M describes the spatial and long-time modulation of the NlW-field. Inserting this 
form in (|A 5|) gives the leading order depth 

r] = — iM^e~ lt + c.c., (A 8 ) 

where c.c. denotes complex conjugate. The evolution equation for M is found at the 
next order by eliminating resonant terms and the dimensional version of the resulting 
equation is then m- Eq. m is also found for continuously stratified flows when 
the geostrophic flow is barotropic. In this case, it applies to the projection of the NIW 
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amplitude onto a single vertical mode, with g'H/f replaced by fr\, where rd is the 
deformation radius of the vertical mode. Mo re details on the deriva tion (in the continuous 
stratification case) can be fou nd in Young fc Ben Jelloull (Il997l). Note that Eq. (12.ip 
differs from that obtained by iFalkovich et al. | (jl994h and iReznik et al\ (|200lh for the 
shallow-water model. This is because they consider a single-layer model in which the 
geostrophic flow is balanced by a sloping free surface. Our assumption of barotropic 
geostrophic flow and consequent absence of interface slope is more relevant to the oceanic 
context where geostrophic flows typically have vertical scales much larger than the mixed- 
layer depth. 


Appendix B. NIW energy in the absence of a background flow 

The energy associated with the linearized reduced-gravity shallow-water system <EU- 
(TOli in the absence of a flow {U = V = 0) is 

E = JJ ^(u 2 + v 2 + erj 2 )dxdy, (B 1) 

using the non-dimensionalisation of Appendix [Aj As expected, (jB ip indicates that NIWs 
have much more kinetic than potential energy. Inserting expansion 422J and solutions 
(lA7l) and 4M| into (ED gives E = Eq + cEi -f- 0(e 2 ), where 

E 0 = fE\M\ 2 dxdy and E x = if (2|M ? | 2 + Me~ u Ul - M%e~ 2it + c.c.^jdxdy. 


Eq is the NlW-kinetic ene rgy appearing in (12.51) . Because U\ varies as e lt (since secular 


terms were removed, see I Young fc Ben Jelloull ( 19971 ). equation (2.25)), the fast-time 


average of E\ is just the first term, which is clearly proportional to I 2 . 

Thus, for i/j = 0, conservation of total energy averaged over fast-time gives (ESI at 
leading order and (|2.6p at the next order. Note that (|2.5p ~ (|2.6p are exact conservation 
laws for the YBJ model but only adiabati c invariants, i.e. a pproximate conservation laws, 


for the parent shallow-water model (e.g., Cotter fc Reich 2004) 
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